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Recent de Haas-van Alphen (dHvA) experiments on high-Tc compounds have been interpreted 
using Lifshitz-Kosevich (LK) theory, which ignores many-body effects. However in quasi-2d systems, 
interactions plus Landau level quantization give strong singularities in the self-energy E and the 
thermodynamic potential Q. These are rapidly suppressed as one increases the c-axis tunneling 
amplitude t± and/or impurity scattering. We show that 2d-3d crossover and interaction effects 
should show up in these experiments, and that they can lead to strong deviations from LK behaviour. 
Moreover, dHvA experiments in quasi-2d systems should clearly distinguish between Fermi liquid 
and non-Fermi liquid states, for sufficiently weak impurity scattering. 
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By tradition de Haas-van Alphen (dHvA) experiments 
are interpreted using Lifshitz-Kosevich (LK) theory, in 
which magnetization oscillations probe directly the quasi- 
particles at the Fermi surface (so that in a non- Fermi liq- 
uid (NFL), with zero quasiparticle weight on this surface, 
LK theory implies no dHvA oscillations at all). Where 
applicable, LK theory allows unambiguous measurement 
of Fermi surface cross-sectional areas, Fermi surface scat- 
tering rates, and Fermi surface band masses pQ. 

Even in 3d, LK theory is not strictly valid because 
of interactions [H [3] ; these cause "Engelsberg-Simpson" 
(ES) deviations from LK, which are seen in experiments 
[4]. In 2d, the mere existence of the Fractional Quantum 
Hall Liquid (FQHL), even when the interaction strength 
V -C fvjj c , shows that Fermi liquid (FL) theory must 
break down in a field, provided impurity scattering is 
weak [5] (ie., once lu c t 1, where oj c is the cyclotron 
frequency and r an impurity scattering time). 

Thus the dHvA experiments recently performed in 
high-Tc systems [H] create a clear paradox. Impurity 
scattering is weak (it must be for a dHvA signal to be 
seen) and the c-axis tunneling amplitude t± is very small 
(in YBa2Cu307_5, t± ~ 15 K is found for 6 = 0.5): thus 
hui c > t± and the system is reaching the 2d limit. And 
yet it is claimed that the data can be fit using LK theory 
[6 . Similar LK analyses have been made for other quasi- 
2d systems [Tj [5]. Since LK theory must break down 
for genuinely 2d systems if oj c t 3> 1 and correlations are 
strong, this raises several important questions: 

(a) How can one generalise dHvA theory to include 
interactions in quasi-2d systems; and how should dHvA 
data then be analysed? 

(b) What kind of oscillations will be shown by NFL 
systems; and can one tell the difference between FL and 
NFL states from dHvA experiments? 

To address these questions, we first analyze the 1- 
particle Green function Q and the thermodynamic po- 
tential 17 for a quasi-2d system, with t±/Hcu c assumed 
arbitrary (but t± <C fi, the chemical potential). When 



interactions are added, we find highly singular behaviour 
in Q. When t± <C h\o c and uj c t 3> 1, these singularities 
imply a complete breakdown of standard Fermi liquid 
theory. However we still find dHvA oscillations, although 
not of LK form. To illustrate these results we compute Q 
for 2 examples; a NFL with singular forward scattering 
interactions, and a FL of band electrons interacting with 
nearly antiferromagnetic spin fluctuations. We find clear 
3d-2d crossover effects as Huj c exceeds t± , and departures 
from LK behaviour whose form depends strongly on the 
nature of the many-body interactions. Neither LK the- 
ory, nor its "ES" generalisation 3 , apply strictly unless 
huj c < tj_ and/or lo c t <C 1; neither condition is satisfied 
in experiments. We find that dHvA experiments ought 
to be able to distinguish FL from NFL states. 

(i) Singularities of Q: The form of the dHvA oscil- 
lations can be found from either the spectral function 
TmQ(e), or directly from fi. In 2d, the Landau levels are 
massively degenerate, and ImQ v (e) oc S(e — e„) where 
e u is the v-th Landau level energy; interactions desta- 
bilize this degeneracy, and so have a singular effect on 
G(e). However any impurity scattering or c-axis tunnel- 
ing tends to suppress this singularity. Although the an- 
alytic structures of Q(e) and fl are now understood for 
neutral 2d fermions [9 in a field (ie. without Landau 
quantization) , there are no general results when one has 
both Landau quantization and interactions |10j . How- 
ever, we can derive results for particular models. Here 
we discuss 2 simple models involving quasi-2d band elec- 
trons, with dispersion = e{k x ,k y ) — 2t±cos(k z a) — fx, 
where t± -C These couple to low-energy fluctuations; 
in a finite field, the lowest-order "1-fluctuation" graph 
for the self-energy takes the form 

£„(fc z ,z)=^^ / — |A^(q)| 2 Xm X (q,a;) (1) 
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where x(q, uS) is the fluctuation propagator, f u — f(e vqz ) 
is the Fermi function for electrons in the ^-th Landau 
level, n(u) the Bose function, and the matrix element 
Aj/j/' (q) > between Landau states v, v 1 and the fluctua- 
tions, incorporates the fermion-fluctuation coupling g q . 
When (i > uj c , |A^^/(q) | 2 ~ g%(m/2fi,)iw c /Trq. 

At this time there is no consensus on a model for high- 
T c superconductors (indeed the central issue is whether 
they are FL or NFL); and other strongly-correlated quasi- 
2d systems are quite complex. Thus, instead of present- 
ing numerical calculations for a specific experimental sys- 
tem, we address the general questions posed in the intro- 
duction by analysing two widely studied models of strong 
correlations in quasi-2d systems: in zero field these de- 
scribe a FL and NFL respectively. 

We begin by discussing the self-energy, which for a 
quasi-2d system can be written near the Fermi surface as 
E(z) = E(z) + E osc (z), where E(z) is non-oscillatory in 
1/B, and the oscillatory part 

E osc (z) = 2£(-l) r E r (z)Jo (^ r J^) sin 

r— 1 ^ c ' ^ 

(2) 

The Bessel function J in this expression comes from in- 
tegrating over q z . 

Model (a) Spin fluctuation model: This well-known 
model [11] has 2d lattice fermions with dispersion 

e(k x , k y ) — —2£o(cos k x + cos k y ) — At\ cos k x cos k y (3) 

and coupling t± between planes; the fermions couple to 
antiferromagnetic spin fluctuations, with propagator 



x(q,w) 
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via a coupling g q = g. The wave- vectors Q = (±7r, ±7r). 
In zero field this model, with or without vertex cor- 
rections |12j . gives FL behaviour, with a Green func- 
tion having finite residue z^ F (/x) at the Fermi surface, 
and a self-energy E(w) with a 2d FL form (ie., with 
lZeE(ui) = (1 — m/m*)uj and ImH(w) oc lu' 2 (1 + In a;)). 

In a finite field, E(w) can be evaluated analytically, but 
the expression is extremely lengthy |13j . The essential re- 
sult is shown in Fig. [T] Landau quantization introduces 
a "step-like" behaviour in <9ZmE/de, with corresponding 
singularities in 7£eE(e), at e = e„. Notice how rapidly 
this singular behaviour is suppressed by interplane hop- 
ping - it is almost invisible once t± ~ hto c . Impurity 
scattering has a similar effect (not shown in Fig. [T]). 

Model (b) Non-Fermi liquid model: We now couple the 
band electrons to fluctuations with propagator [H] : 



X<1 



(5) 
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FIG. 1: The imaginary part ImS(f) of the self-energy, as 
a function of t±/hui c , at T = 0. (a) shows ImdT,/de for 
the spin fluctuation model; we fix g q — g = 0.58 eV, xo — 80 
states/eV, £ = 2.5a and usf = 10 meV. (b) shows ImS (with 
no derivative) for the non-Fermi liquid model, assuming s — 3; 
we fix K s=3 = 0.013(fitJ c ) 2/3 ^ 1/3 with fx = 6000 K. (c) shows 
the effect of impurity scattering on the NFL model; we plot 
ImS(t) for different values of u> c t, assuming t±/fkJ c ~ 0.75. 



field self-energy has the NFL forms E(e) ~ elne (for 
s = 2) and S(e) ~ (if2o/e)se (for s = 3), so that 2k(e) — * 
on the Fermi surface. In a finite field, E„(fc z ,z) can 
again be evaluated analytically in the form ([2]), with the 
coefficients E r (z) taking the interesting T = form : 



E r (z) 
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where s is a dynamic scaling exponent, with 2 < s < 3, 
using a fermion fluctuation coupling g q — K s . The zero 



(and a more complicated finite T form), where S2 is a 
Lommcl function |15j . and Z r — 2ixrz /frw c . Now the 
singular behaviour in E is far more pronounced; again, 
it is eliminated by switching on tj_ (Fig. |TJb)) , or by 
impurity scattering (calculated in Fig. [ljc) in a self- 
consistent Born approximation). 

We see that both models show singular behaviour of 
E(e) as a function of e, implying similar behaviour for the 



3 



quasiparticle weight Zfc(e). At the Fermi energy, E(e = fi) 
will then show the same singular behaviour as a function 
of B, periodic in \/B. Strictly speaking, this means a 
breakdown of FL theory for both models, but much more 
strongly for the NFL system. Because these singularities 
are rapidly suppressed by both inter-plane hopping and 
impurity scattering, this breakdown will only be clearly 
visible when t±, h/r <C Hu) c . 

(ii) Thermodynamic potential fi: If "crossed 
graphs" can be ignored in we can write an expres- 

sion for f2 in terms of Q [TBI : 



n 



i, 

0' 



-Tr In 



(G + Gosc) 
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where Q is the non-oscillatory part of Q. This expres- 
sion resembles the classic Luttinger/ES expression [5] 
for fl, except that the latter drops Q osc from This 
is justified in 3d, but not in 2d |16j ; in the quasi-2d 
case it is only justified if tj_ ^> huj c . From we find 
fl = £l + 2Y l 7 =1 (-l) r n r cos(rhA F /eB), where ft is the 
non-oscillatory part of fi, and 



— V [ J 2 ( 47rrt± 

~h 2 i3 ^ 1 ° I hw c 

H n>0 v c 



Cr(^„) 



(8) 



hjjj c ( 47trt± \ 2nr 

7> — J o ~Z CXp(- — (W ri 

27rr \ nu; c / n,w c 



CK))] 



where the Cr(w„) = iE r (za; n ) are real and positive. Equa- 
tion (|8| reduces to the Luttinger/ES expression for Q if 
we drop the first term, and if in the second term we use 
only the non-oscillatory part of C( w n) = iE{io} n ). 

It further reduces to LK if we assume — * (1 — 
m/m*)cu + i/2r, ie., a mass renormalisation and scat- 
tering rate both independent of energy. Clearly (Fig. [TJ) 
the oscillatory part of £ must not in general be neglected. 

(iii) Oscillatory Magnetisation: We write the 
magnetisation at constant chemical potential M^(B) — 
- dfl/dBl^ in the form M = M(B) + 2 £ r (-l) r M r , 
where M is the non-oscillatory part (M(B) at constant 
N is found by making a Legendre transform |17lL Dif- 
ferentiating pj, we get M r = M{ + MJ, with 
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The key point here is that if £ contains strong oscilla- 
tions with energy, these translate into a very strong new 
oscillatory contribution to M osc . 

Equation (JoJ) yields a very rich variety of forms for 
M{B), depending on the two parameters t±/Hu c , uj c t, 
and on the form and strength of the interactions. We 
have no space here to discuss the whole parameter range, 
but we can summarize the key features: 




temperature (in K) 

FIG. 2: lnM r for r = 1 against T (the "mass plot"), for the 
spin fluctuation model: (a) Plot for different fields, assum- 
ing t± — 15K but no impurity scattering (when B — 55 T, 
hu c = 20 K); (b) Plot for different impurity scattering rates, 
assuming t±/hu} c = 0.6. The dashed lines are fits to LK. 



1. Clear departures are seen from LK theory, even in 
mass plots (Fig. [2]), for NFL systems and even for FL 
unless the fluctuation energy scale u>sf ^ ^w c , and/or 
lu c t < 1 . Without interactions, only the 2nd term (oc Jo) 
survives in ([8]) ; this term is well-known in LK theory |18j . 
With interactions, the two terms compete - as the field- 
induced singularities in S(e) (and hence in C( w n) become 
stronger, the 1st term in fl r (oc J5 2 ) increases, and for 
NFL it can dominate the 2nd term. The much stronger 
singular structure in S(e) means that NFL have much 
stronger departures from LK than FL. 

2. The form of M(B) depends strongly on t±/Huj c . 
This gives a remarkable structure in field plots (Fig. |3|, 
which is eliminated by strong impurity scattering (Fig. 
|2jb)), or by removing strong correlation effects. 

3. Short-range impurity scattering strongly suppresses 
the singular structure from interactions once lo c t < 1 
(see Fig. [2^b)). However, curiously, it affects M[ and 
M£ rather differently; M{ decreases exponentially with 
1 jio c T (a la Dingle) but M£ decreases approximately as 
a power law. More refined analysis of the effect of scat- 
tering off impurities and off small angle scatterers (like 
dislocations) is certainly necessary for this problem. 

We see that interactions have profound effects on the 
quasiparticles and the thermodynamics of conducting 
systems in high fields, for quasi-2d systems. These ef- 
fects are rapidly removed by interplane coupling (once 
tj_ > huj c ) , and even more rapidly by impurity scattering 
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FIG. 3: The dHvA component M r for r = 1, against 1/B 
(again, if B = 55 T, h~uj c = 20 K), for different values of 
t± (measured in K), with no impurity scattering; we assume 
T = 1 K. (a) for the spin fluctuation model; and (b) for the 
non-Fermi liquid model with s = 3 (the topmost curve in (b) 
has been rescaled by a factor 0.5). 



(once uj c t < 1). The models we have used are of course 
rather simple (although very widely used in the litera- 
ture); but our main results are not crucially changed by 
eg., adding vertex corrections. 

Consider now the experimental situation. Experiments 
on YBCO fall precisely in the crossover regime, with 
t± ~ 15 K, and 15 K < hui c < 30 K. It is not yet pos- 
sible to compare the experimental fits [6 J on YBCO and 
Tl-2201 with the theory here, because these fits have not 
included the Jo term (which already exists in LK the- 
ory [IE])- It b e extremely interesting to have fits to 
different strong-correlation models - and to discriminate 
between FL and NFL models. We note that absence of 
the Jq term in Q would indicate the underlying state is 
FL (but NFL if the Jq term is strong) . It will also be in- 
teresting to look more closely at other strongly-correlated 
quasi-2d systems in high fields - where few departures 
from LK have been found so far. Finally, note that any 
experiments sensitive to the singular structure we find in 
Q should show interesting effects. Obvious examples are 
c-axis tunneling and SdH experiments in very high fields, 
but a generalisation of the foregoing to a transport theory 
will be required. 
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